Mon. Not. R. Astron. Soc. OOP. [THTSl (2013') Printed 4 April 2013 (MN WI^ style file v2.2) 
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ABSTRACT 

We present a weak lensing analysis of the cluster of galaxies RXC J2248. 7-4431, a 
massive system at z = 0.3475 with promine nt strong lensing features covered by 
the HST/CLASH survey (iPostman et al.ll2012[ ). Based on UBVRIZ imaging from the 
WFI camera at the MPG/ESO-2.2m telescope, we measure photometric redshifts and 
shapes of background galaxies. The cluster is detected as a r nass peak at 5a si, 
nificance. Its density can be parametrised as an NFW profile ([Navarro et al.l [199' 
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1 INTRODUCTION 

Clusters of galaxies, such as RXC J2248.7-4431 (cf. Figure[T|) 
studied in this work, are the most massive gravitationally 
bound objects that have formed in the universe until today. 
What makes them interesting for cosmology is that they 
lie at an intersection of the two potentially most important 
unresolved questions: their mass content is dominated by 
Dark Matter - and their formation and evolution is strongly 
influenced by the interplay of matter density and Dark En- 
ergy. For this reason, studying clusters of galaxies is also a 
powerful probe of cosmological parameters an d models (e.g. 
Wang fc Steinhardd 1 19981 : [Holder et al] 1200 ll : IWeller et al.l 



20021 ) 



Perhaps the most important property of a cluster of 
galaxies with respect to cosmology is its virial mass, which 
can be determined in several different ways. The num- 
ber and the velocity dispersion of cluster member galax- 
ies are related to the total mass of the system. One 
can use larger fractions of the overall mass as a proxy 
by observing the hot intra-cluster gas, which emits X- 
ray radiation (e.g. [ B5hringer et al.|[2004l : ISahlen et al.ll2009l : 
iPiffaretti et al.l I2OIII ) and changes the cosmic microwave 
background spectrum due to invers e Compton scattering 
(Sunyaev-Zel'dovich (SZ) effect, cf. ISunvaev fc Zel'dovicbl 
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Figure 1. BRI colour image of the central 200x200 arcsec of 
RXC J2248. 7-4431. Clearly visible are the brightest cluster galaxy 
and several gravitational arcs, along writh a high concentration of 
yellowf cluster member galaxies. 
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1 19721 : iHaiman et allbOOll : iBattve fc Weiieill2003l ). The lat- 
ter observables can be relat ed to mass by astrophysical 
modelling or self-calibration (|Hull2003l : iMaiumdar fc Mohn 
|2004 ). 

Despite this, however, weak gravitational lensing is a 
valuable ingredient since the tangential alignment of back- 
ground galaxy images is directly proportional to overdensity 
of all cluster matter - luminous or dark - alike and is insen- 
sitive to the astrophysical state of the cluster. This allows 
for unbiased mass measurements of single clusters and im- 
proved calibration of other mass-observable relations (e.i 
lAUenet al.ll2002l : lOkabe et al.ll2010l : iHoekstra etal] 12011 
Furthermore, gravitational lensing straightforwardly allows 
the probing of additional properties of the density profiles of 
dark matter haloes, such as their concentrations, for which 
predictions in a cosmological model can be made and from 
which additional constraints can be drawn (see for instance 
the Cluster Lensing And Supernova survey with Hubble 
(CLASH), in which RXC J2248. 7-4431 is also observed, cf. 
IPostman et al.ll2012h . 

In this work we analyse the weak lensing effect of 
the cluster RXC J2248. 7-4431 based on background galaxy 
shapes and photometric redshifts measured from UBVRIZ 
multi-band imaging by the Wide-Field Imager (WFI) on the 
2.2m MPG/ESO telescope at La Silla. In Section 2, we in- 
troduce the basic properties of the data used and our data 
reduction. Section 3 summarises previous optical. X-ray and 
SZ observations of RXC J2248. 7-4431. We give an overview 
of the photometric analysis, including photometric redshifts, 
cluster member photometry and morphology in Section 4. 
Technical aspects of our weak lensing measurements are dis- 
cussed in Section 5. The analysis of the weak lensing effect 
of RXC J2248.7-4431 is presented in Section 6. Results for 
a second cluster found in the field of view at z « 0.6 are 
shown in Section 7. We summarise our results in Section 8. 

All numerical values given in this work are calculated 
for cosmological parameters Ho — 72km s~^ Mpc"^ and 
Q,m = 1 — Ha = 0.27. Where applicable, measurements 
from the literature have been converted to this cosmology 
as well. We denote the radii of spheres around the cluster 
centre with fixed overdensity as rAm and tac, where A is 
the overdensity factor of the sphere with respect to the mean 
matter density pm or critical density pc of the universe at 
the cluster redshift. Masses inside these spheres are labelled 
and defined correspondingly as MAm = A • ^r\^ ■ pm and 
Mac = A • f ri, • Pc 



2 OBSERVATIONS AND DATA REDUCTION 

This analysis is based on observations made with the Wide- 
Field Imager (WFI) on the 2.2m MPG/ESO telescope at 
La Silla. The sensitivity of the instrument and the filters 
available spans all optical wavelengths and with its field of 
view of 33 x 33 arcmin it is well suited for weak lensing cluster 
analyses with photometric redshifts. 

The observations in U (#877), B (#842), V (#843), 
R (#844), I (#879) and Z (#846) bancfl used in this work 
were taken in the years 2009-2010. Details of integration 
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Figure 2. Cumulative exposure time as a function of limiting 
seeing in WFI frames of RXC J2248. 7-4431 used in this work. 
Bold lines show the bands used for shape measurement, where we 
have applied seeing full-width at half maximum (FWHM) cuts at 
0.9 arcsec (1.0 arcsec) in R (V and I) band, respectively. 



time as a function of limiting point-spread function (PSF) 
full-width at half maximum (FWHM) are shown in Figure[2l 
The excellent depth and seeing in R band (with over 9 hours 
of exposure time at sub-arcsecond seeing, yielding a coadded 
image with 0.8 arcsec PSF FWHM and limiting magnitude 
of rn,i{_voga,iim = 26.7 for a 5a detection inside a 1 arcsec 
diameter aperture) make it our primary lensing band, but 
also in V (0.9 arcsec seeing with m\/_vega,iim = 26.4) and I 
band (0.9 arcsec PSF FWHM with mj,vcga,iim = 24.4) shape 
measurements of background galaxies are feasible. 

We perform de-biasing, fiat fielding, masking of bad 
pixels in all bands and fringe pattern corr ection in the I 
and Z band using the Astro- WISEJj pipeline (IValentiin et al.l 
|2007|). Background subtraction, final astrometry and co- 
addition of su itable frames i s done externally us ing SCAMIjj 
(|Bertinll2006l ) and SWarjQ iJBertin et al.ll2002h . The central 
part of a colour image based on these co-added frames is 
shown in Figure [l] 

For photometry, we take observations of standard star 
fields in B, V and R to fit zeropoints for each CCD individ- 
ually and extinction coefficients globally. From nights which 
we confirm to be photometric according to their standard 
star measured zeropoints we build a photometric B, V and 
R stack. We build deeper stacks including frames from all 
available nights with relative photometry fitted by SCAMP 
and match the magnitudes measured from these stacks to 
the photometric ones in order to find zeropoints for the deep 
stacks. Alternatively, fixing the R band zeropoint, we find U, 



^ http://www.astro-wise.org/ 
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B, V, I and Z zeropoints by minimising residuals in colour- 
colour diagrams with respect to the stellar library of lPicklea 
(| 19981 ). The stellar locus zeropoints found in this way for 
B and V agree within 0.02 mag with the photometric ze- 
ropoints found by matching the two versions of the stacks, 
confirming that our photometric pipeline provides consistent 
results. 



3 PREVIOUS WORK ON RXC J2248.7-4431 

The cluster RXC J2248.7-4431 studied in this work is also 
known as AbeU S1063 or MACS 2248-4431. Here we give an 
overview of all literature on the cluster published previously 
to this work, including detections and redshift estimates, X- 
ray observations an d SZ mea surement. We discuss the more 
detailed findings of I Gomez et al. ( 20121) and compare them 

to our own analyses in Section [71 

RXC J2248. 7-4431 was first listed bv lAbell et all (|l989l ') 
with a background-corrected galaxy count of 74. It was 
independently detecte d by the ROSAT All-Sky Survey 
l|de Grandi et al.l 19991 ). who quote a redshift estimate of 2: = 
0.2 52, which was, how ever, only weak ly determined based on 
the lAbell et all (|l989l ') distance class. ICruddace et al.l (|2002l ') 
give its redshift as z = 0.1495 bas ed on private communi - 
cations with H. Andernach (1989). iBohringer et al.l (|2004l ') 
finally quote a spectroscopic redshift of Zd = 0.3475, which 
is what we adopt for this work. This is confirmed by spec- 
troscopy of 51 cl uster mem bers with a mean redshift of 
(iGomez et al. 2012). 



Z — U.J^Dl_nnnii 



iGomez et al.l ( 
persion of to be (t„ 
a mass of M200C 



20121 ') find the spectroscopic velocity dis- 



1660ti5okm s ^ , which corresponds to 
421 q^ • 10^*Mq according to the relation 
of lEvrard ct al." (2008'). 

iCruddace et al. (2002) quote an X-ray temperature 
of 7.823 keV f r om th e ROSAT All-Sky Survey (RASS). 
iMaughan et al.l ([2003) give the X-ray temperature within 
7?5oo as II.iIq gkeV and parametrise the cluster profile as 
slightly elliptical with 1 — b/a = 0.2 based on Chandra 
ACIS-I data. Both values are based, howev er, on earlier er- 
roneo us cluster redshifts (2 = 0.1495 for Cruddace et al.l 
20021 and 



0.252 fo r iMaughan erall 120081 ). The 



Planck Collaboration et al.l (|2011bl ) quote an X-ray mass 



based on their follow-up with XMM-Newton of Msooc = 
(12.25 ±0.21)-10^'^MfT) w hich they use to calibrate their SZ 
MOE. IComis et al.l l|201ll ) determine M2500C = (5.3 ± 2.6) ■ 
1O"M0 from Chandra X-ray data and lGomez et all (|2012l ) 
give a consistent value from independent data reduction of 
the same data of M2500C = (6.0 ± 1.6) • 10^" Mq. 

IPlagge et al.l (|20ld ) parametrise the SZ profile mea- 
sured with the South Pole Telescope (SPT) with a /3 pa- 
rameter of (0.86 ± 0.02) arcmin at a scaling of AT ~ 1 mK. 
IWillia mson et al.l ||201H ) quote an SZ signal-to-noise ratio 
(S/N) of 20.7 and a SZ mass of M2oo^ = (28.2± stat 3.6±sys 
9.3) • 1O"M0. The lPlanck Collaboration et al.l (|2011al ) de- 
tect the SZ effect of RXC J2248. 7-4431 at 13.93cr signifi- 
cance. From the Planck SZ observ able and scaling relation 
(jPlanck Collaboration et al.ll2011a ) the SZ mass is M500C ~ 

(11. 5 ±,tat 2.6 ±sys 0.5) ■ VS'^Mp 
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Figure 3. Three-dimensional galaxy density as a function of pho- 
tometric redshift in the WFI field (black, solid lines) and inside a 
circle of radius 100 arcsec around the brightest cluster galaxy of 
RXC J2248. 7 -4431 (red, dashed line s). The cluster at a redshift of 
2ci = 0.3475 IjBohringer et al.ll2004r can be clearly seen as a peak 
in this smaller field. The green, dotted line shows objects within 
a radius of 100 arcsec around the BCG of the z ft: 0.6 cluster 
described in Section |8] 



iGuzzo et all (|2009i ) and lPlagge et all (|2010D show X-ray 
and SZ imaging of the cluster, respectively. The cluster is 



also covered by th e CLASH project (jPostman et al.ll2012l : 
iMonna et aDl2013l ). 



4 PHOTOMETRIC ANALYSIS 

4.1 Photometric Redshifts 

The multicolour catalogue creation and photometric 
redshift estimat i on fo llows the procedure described in 
iBrimioulle et al.l (|2013l ) . Here we only give a brief overview. 
We convolve all data with a Gaussian kernel to match the 
seeing to the band where it is largest (in this case the U 
band). This equalisation of the PSF is required for reli- 
able aperture colours in the different filters. We then run 
SExtractoijJ in dual-image mode to extract fluxes, includ- 
ing weight images and masks of bad areas in the detection 
and extraction frame. We detect the objects on the uncon- 
volved R-band with a S/N threshold of 2o" on at least four 
contiguous pixels. Flux and magnitude information is ex- 
tracted from the convolved images. We correct for extinc- 
tion and zeropoint accuracies by comparing the stell ar locus 
in col our-colour-diagrams to the stellar library of IPickled 
(|l998l ). We then use the pho tometric template-fitting algo- 
rithm of iBender et al.l l|200ll ) to estimate photometric red- 
shifts. 

The density of galaxies in redshift space is shown in 
Figure O The cluster is clearly visible as a redshift density 
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Figure 4. Density of galaxies in our piiotometric redshift cata- 
logue in a map centred on the BCG position, averaged in angular 
boxes of 4x4 arcmin size. Contour lines show density per arcmin^ 
per redshift interval for sources with [z — z^i! < 0.04- {1 + z^i). The 
background colour shows the overall density of sources above the 
redshift cut of Eqn. [T] The depletion of background galaxies at 
the cluster position is an indication that our background sample 
is pure. The high background density observed in the lower left 
of the map is due to the system discussed in Section |8] 



peak around the spectroscopic Zd = 0.3475, which is also 
true for the second cluster at a mean photometric redshift 
of z ~ 0.6, described in detail in Section [8] 

In our weak lensing analysis we only consider back- 
ground galaxies with a minimum photometric redshift of 



> 1.1 -Zci-f 0.15 ^0.53 



(1) 



We include objects up to a maximum photometric redshift of 
z — 4. We verify that this excludes redshift mismeasurement 
of cluster galaxies into the background sample by means of 
redshift density maps at the cluster redshift and for objects 
selected by the background cut of Eqn.[l]in FigurelD Indeed, 
the cluster is seen as a highly significant overdensity of ob- 
jects around z^i, while the density of background sources 
near the cluster is smaller than in the field. This is expected 
owing to the area covered by cluster hght and indicates that 
background contamination is not a significant issue, as ex- 
pected for the red cluster galaxies whose redshift can be 
determined photometrically to good accuracy. 



4.2 Cluster member SEDs 

We investigate photometric properties and spectral energy 
distribution types (SEDs) of the cluster member galaxies of 
RXC J2248. 7-4431 and the second cluster at z Ri 0.6 (cf. 
Section [8|. 



' " x' J'T 

X. X 







Figure 5. Colour-magnitude diagram of galaxies within 5 arcmin 
projected separation from the BCG, using B-V aperture colours 
and R band MAG_AUTQ magnitudes and showing all objects (black, 
small crosses), objects with red spectral energ y distribution types 
(cf. Dahlon ot al. 20051: iBrimioulIe et al.ll2013h (red, large crosses) 
and marking objects in a redshift slice \z — z^l ^ 0.04 ■ (1 -|- z^i) 
around the cluster by a green open square. 



The brightest cluster galaxy (BCG) of RXC J2248.7- 
4431 has a photometric redshift of 0.4 and an absolute mag- 
nitude of Mfl^vcga = —24.81. Its counterpart in the z ~ 0.6 
cluster has a photometric redshift of 0.66 and an absolute 
magnitude of Afij^voga = —25.1. 

Cluster members galaxies are preferentially of early 
type and, at a given redshift, lie on a relatively tight se- 
quence in colour-magnitude space, the red sequence. We 
plot B-V aperture colours, which are particularly indicative 
of the D4000 break, against R band MAG_AUTD magnitudes 
in Figure O The red sequence is clearly detected and dom- 
inant among cluster members. We measure the fraction of 
red galaxies on a sample of cluster members selected by their 
photometric redshift |z — Zci| ^ 0.04 • (1 -|- z^i) for both RXC 
J2248. 7-4431 and the second cluster found at z ^ 0.6. Fig- 
ure [S] shows the fraction of these galaxies wh ich are classi- 
ficd photometr ically as havin g a red SED (cf. iDahlen et al.l 
2005; Brimiou Ue et al.ll2013r ). Nred/Ntotai, as a function of 
projected separation from the cluster in units of its r2oom 
according to our best-fit NFW models (see Sections 16.51 and 
[8| and as a function of depth. 

The latter is fixed relative to the characteristic lumi- 
nosity of the Schechter function, L*. L* , or in this case m* , 
is estimated by fitting a Schechter function to the statisti- 
cal ly background corre cted cluster luminosity function (as 
in IZenteno et al.l I2OIII ) . The limiting depth for the fit is 
choosen at m* + 4, where m* is found iteratively. The off- 
set of 4 limits the magnitude range so it is well sampled 
for both clusters. In this way we find characteristic mag- 
nitudes of m|j,vega = 19-3to;g ('^fl.voga = 22.2lj 7 for the 
z Ri 0.6 cluster) for apparent R band and M* ab ~ — 23.4^0 3 
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-23.9 



+0.9N 

i.oJ 



for absolute magnitudes in SDSS r' 



band. 

We conclude that approximately 50 per cent of cluster 
members are red galaxies near the core of RXC J2248.7- 
4431, and that this fraction continuously drops towards the 
outskirts, where it is below 20 per cent. For the higher red- 
shift system, the fraction of red galaxies in the cluster is 
signific antly lower (which is kn own as the Butcher-Oemler 
Effect. iButcher fc Oemieil|l97 8l). The fraction of red galax- 
ies decreases both towards fainter magnitudes, larger sepa- 
rations from the core and higher redshift. 



4.3 Mass from richness and luminosity 

We can find an independent estimate of cluster mass by 
means of mass-observable relations of richness and abso- 
lute lumin o sity. To this end, we apply the prescription of 
IWen et al.l (|2010l ). For our cluster member catalogue we se- 
lect luminous galaxies inside a projected radius of IMpc 
around the BCG within \z — Zci\ ^ 0.06 • (1 -I- Zd) and with 
magnitudes M^' ^ — 21, where we u se absolute magnitudes 
in SDSS r' band (|Gunn et al.lll998l 'l. The size and r' band 
luminosity of the sample is 64 and L^' = 3.0 • 10^^ Lq. Af- 
ter subtracting the mean number and luminosity of similar 
object selected from a region of the field more than 3Mpc 
from the BCG, we find a background corrected richness of 



i? = 54 and luminosity of L^- = 2.6 ■ lO^^L, 



^ 



Applying the scaling relations of lWen et al.l (l2010|), this 
yields 



logMioic/(lO'*/i-'M0) 



= (-1.57±0.12) + 1.55 log i? 
= 1.12 ±,tat 0.09 ±sys 0.12 (2) 

from richness as the mass proxy with systematic and statisti- 
cal uncertainty. Alternatively, one can use the net luminosity 
to obtain 

logMioic/(10^*/i~^MQ) = (-2.03 ±0.06) 

+ 1.49 log L,7(10'''/i"^Lo) 
= 1.14±,tat0.09±sys0.06.(3) 

Since the errors of these two estimates are highly correlated 
we use the less uncertain Eqn. [S] which yields a mass of 



Ml. 



19 ± 6 ■ W'*M, 



0- 



5 WEAK LENSING MEASUREMENT 

Weak gravitational lensing changes the ellipticities of the 
images of galaxies in the background of massive structures. 
Its measurement therefore requires the determination of pre- 
seeing galaxy shapes. In this section we describe our shape 
measurement, while the analysis of the signal is presented 
in Section [B] 



5.1 Model of the Point-Spread Function 

Any estimate of the pre-seeing shape of a galaxy requires 
knowledge of the point-spread function (PSF). For a suc- 
cessful and unbiased weak lensing analysis it is most cru- 
cial to model the point-spread function accurately (cf. 
iKitching et al.ll2013l ). since mismatches in ellipticity and size 




2 4 6 

Flux radius [pixel] 

Figure 7. Selected stars for PSF modelling (red) and sources 
initially accepted for shape measurement (green). 
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Figure 8. Autocorrelation of PSF model ellipticity residuals (-Di) 
and cross-co rrelati o n of residuals and measured ellipticity (-D2) 
according to iRowj 1 I2OIOI ) for a 5th order polynomial PSF model 
in the R band co-add. While lower orders show signs of underfit- 
ting (Di and/or D2 significantly non-zero on some scales), higher 
orders do not improve the model further and show signs of over- 
fitting (D2 significantly negative on some scales). 



© 2013 RAS, MNRAS OOO.ITlfTsl 



6 D. Gruen et al. 




0.8 t 



iO.6 - 



2 0.4 



0.2 - 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r- 




■l-K. 



-K.: 



-f— i-. 



Figure 6. Fraction of galaxies having a red SED type in the rodshift shce and projected vicinity of RXC J2248. 7-4431 (sohd hnes) and 
the second cluster at z ^ 0.6 (dashed lines). The left panel shows the red fraction as a function of radius from the BCG in units of 
''200m of the respective system writh magnitude limits at mj, + 1 (red), m^ + 2 (green), rrip + 3 (blue) and our detection limit (black). 
The right panel shows the red fraction inside r200m of the respective system as a function of magnitude limit, which reaches the overall 
detection limit at m^™ — mjj Ri 4 for the z Ri 0.6 system. Error bars are calculated from the Poissonian noise in the number counts of 
red and blue galaxies in each bin. 



of the PSF model cause additive and multiplicative system- 
atic errors in the sha pe estimate, respectively. 

We use PSFEjfl (|Bertinll201ll ') for determining the pro- 
file of the PSF as a polynomial function of the position in 
the co-added focal plane. To this end, we perform a pre- 
selection of stars according to siz e, S/N and the SExtrac - 
toJJ neural network star classifier (JBertin fc Arnoutal 19961 ). 
The pre-selected stars (cf. Figure [7|) are used to determine 
the PSF model, which is then checked using the following 
diagnostics: 

• mean and scatter of residual ellipticities and sizes to 
check for remaining offsets and quality of the fit, 

• whisker plot of residual ellipticities to exclude failure of 
the fit in specific regions of the focal plane, and 

• two-point autocorrelations of ellipticity residuals (Di) 
and cross-correlation of residuals and measured star ellip- 
ticitie s (D2) to exclude under- and overfitting, as defined in 
iRowd l|2O10h . Eqn. 13-14. 

We verify, in particular, that a star preselection is nec- 
essary for the size of the PSFEx model to match the size 
of the PSF well enough. In R and V a smooth PSF model 
over the whole focal plane is a sufficient description of the 
observed pattern. We find that, similar to other exposures 
from the WFI camera we have analysed in the past, the I 
band PSF is more difficult to model. Only when masking 
the border regions between neighbouring chips and discard- 
ing all stars and galaxies in those regions can we describe 



^ http://www.astromatic.net/software/psfex 

^ http://www.astromatic.iiet/software/sextractor 



the PSF with a simple polynomial dependence on position to 
suffi cient accuracy. For these tests we find that the method 
of Rowel (|201Cr ) is particularly helpful, and use it to deter- 
mine the correct polynomial order of the spatial variation 
(5 (7) for the R (V and I) band, respectively), which match 
the observed PSF without clear signs of over- or underfitting 
(cf. Figure (SI). 



5.2 Shape Measurement 



We 



run an 



of 



(Kaiser et al. 



imp lementation 
1995 (KSB), [Tuppino 



the KSB -1 - me thod 
Kaisej Il997l and 



iHoekstra et al. 1998 ) using the PSFEx PSF model (KSBPS- 
FEx) that has been tested against simulations at intermedi- 
ate to large tangential shears and proven to be viable in the 
cluster shear regime (jYoung et al.l 120131 ). The pipeline in- 
cludes the following preparation steps for the KSB+ shape 
measurement: 

• unsaturated sources with flux radii larger than the stel- 
lar flux radius and zero SExtractor fiags are filtered, 

• postage stamps of 64x64pix size are extracted and 
neighbouring objects according to the SExtractor segmen- 
tation map are masked, 

• the SExtractor photometric background estimate at the 
object position is subtracted from the image in order to com- 
pensate for small-scale background variations insufficiently 
modelled by the data reduction pipeline, 

• bad and masked pixels are i nterpolated using a Gauss - 
Laguerre model of the galaxy (jBernstein fc Jarvia [2003) i 
where we discard objects with more than 20 per cent of 
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postage stamp area or 5 per cent of model flux faUing onto 
bad or masked pixels. 

As the final step of the shape measurement, KSB+ is run on 
the cleaned postage stamp of the galaxy and the sub-pixel 
resolution PSF model at the galax;y position. For details, we 
refer the reader to the papers introducing and extending the 
method. Here we only give a brief summary. 
KSB measure polarisations 



1 



Qii + Q22 



2Qi2 



(4) 



using second moments calculated inside a Gaussian aper- 
ture, 



Qij 



I 



(fei(0)w{\e\)eit 



(5) 



of the surface brightness distribution of the galaxy I{G) with 
a Gaussian weight function w^(|0|) centred on the galaxy cen- 
troid. In our implementation, the weight function for mea- 
suring the galaxy and PSF moments is scaled with the mea- 
sured half-light radius of the observed galaxy. 

Cluster weak lensing analys es require the ensemble 
meas urement of reduced shear (cf. iBartelmann fc Schneideij 
I2OOII . p. 60). In the presence of an elliptical PSF, the linear 
approximation of how observed post-seeing polarisation Go 
reacts to a reduced shear g can be expressed as 



ei + P p + P'g 



(6) 



where e^ is the intrinsic post-seeing ellipticity of the galaxy, 
P^'^ is a 2 X 2 tensor quantifying the response of observed 
shear to PSF polarisation p and P'*' is the shear respon- 
sivity tensor. Inverting {P'^)~^ x ,^ and assuming that 
(P^)~^ei is zero on average because of the random intrinsic 
orientation of galaxies, this yields the ensemble shear esti- 
mate 



(g) = (e) = 



trPT 



Psm 
P 



(7) 



We remove objects with failed KSB-(- measurements and 
objects with trP'' < 0.1 from the final shape catalogue. 
The latter is a requirement owing to the noisiness of mea- 
sured trP'', which o therwise greatly ampl i fies the uncer- 
tainty of Eqn. [7] (cf. Applegate et al. 20121: Newman et al.l 



I2OI3I : ISchrabback et al.ll2007l : iRomano et al.ll2010l . for sim- 
ilar clipping schemes at various levels of trP^). We match 
R, V and I band shape catalogues and take the arithmetic 
mean of all available bands as the individual object shape 
and match against the photometric redshift catalogue to se- 
lect valid background galaxies and assign redshift scalings. 
This and the S/N cut introduced in Section 15.31 leaves us 
with an average of 9 galaxies per square arcminute in the 
background of RXC J2248. 7-4431 with a mean ratio of angu- 
lar diameter distances of (Dds/Da) ~ 0.59 (cf. Section [6. ip . 

5.3 Noise Bias Calibration 

Biases of shape estimators are commonly expressed as a 
multiplicative and an additive term ( Heymans e t al. 2OO3), 
where the additive component is typically due to imper- 
fect correction for the ellipticity of the PSF. The latter ef- 
fect is particularly relevant for cosmic shear analyses, where 
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Figure 9. Multiplicative bias m as a function of signal-to-noise 
ratio (S/N) for our shape pipeline, as measured using image sim- 
ulations with a range of PSF sizes, shears and galaxy properties 
typical for a weak lensing cluster study like ours. The results 
are well-fit by Eqn. [9] with best-fitting parameters A = 0.025, 
B = 0.17 and C = 17 (red line). 



auto-correlations of the shear field are calculated and auto- 
correlations of the PSF ellipticity field would enter the equa- 
tion by means of an additive bias (for a more detailed 
study of the infiuence of shape biases on cosmic shear, see 
lAmara fc Refregieill200a ). All analyses we do in this clus- 
ter weak lensing study, however, deal with shear estimates 
averaged in an annulus, in which case any constant addi- 
tive bias cancels out. Pixel noise, on the contrary, typically 
causes multiplicative biases at the per-cent level over a wide 
range of observational conditions, which wo uld directly en- 
ter our cluster weak lensing model (jMelchior fc Viola 2012 ). 
This and the fact that we do not observe significant additive 
biases on simulated images with elliptical PSF is the reason 
for limiting our calibration to the most important factor, 
namely noise dependent multiplicative bias. 

Shape measurement, primarily because it involves 
division by noisy quantities, suffers from a noise bias, 
which is negligible for bright galaxies but becomes trou- 
blesome for the faint e nd of backgrou nd galaxies in a 
weak l ensing analysis ( Kaisei I2OO0I: iBernstein fc Jarvid 



20021: iMelchior fc Violal |2012| : iRefregier et al.1 I2OI2I : 
Kacprzak et al.l l2012f) . Many studies in the past have 
applied a glob al correction fact o r to their KSB shape 



catalogue (see IClowe et al.l 20121: Schrabback et al.l 120071 : 
iNewman et al.l 120131 : iRagozzine et al.l I2OI2I for a range of 
factors used), but there also have been cases where calibra- 
tion has been matched to galaxy S /N and sometimes also 
size using simulated galaxy images (|Schrabback et al.ll2010l : 
Ivon der Linden et al.ll2012l ). We quantify the magnitude of 
noise bias in our pipeline using image simulations with a 
range of PSF sizes, shears and galaxy properties typical for 
a weak lensing cluster study like ours (|Young et al.ll2013l 'l. 
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If a shear estimator go only has muhipUcative bias m, 
it can be written as 



{l + Tn)-g + N 



(8) 



where g is the true shear and A'^ is a noise term with zero 
mean. We fit m for simulated sets of galaxies and show the 
results in Figure [9] as a func tion of S/N measured according 
to the prescription of (Erben et al|2001i). The multiplicative 
bias is fit well by the functional form 



m = ~A- B- exp(-(S/N)/C) 



(9) 

: 17 



with best-fit values oi A = 0.025, B = 0.17 and C 
down to a signal-to-noise ratio of 10. 

We discard objects outside the r egime this calibrat ion 
was tested on (i.e. where S/N< 10 (JErben et al.l 120011 ) or 
(S/N)iso < 15 as measured inside an isophotal aperture of 
1.5ct significance per pixel over the background). The re- 
mainder of our sample is calibrated with the m of Eqn.|^ by 
multiplying shape estimates by a factor of 1/(1 + m(S/N)). 



6 WEAK LENSING ANALYSIS 

6.1 Introduction 

Cluster weak lensing analysis aims to reconstruct properties 
of the density field of clusters of galaxies from the reduced 
shear they impose on the images o f background galaxies. We 
refer the reader to the review of iBartelmann fc Schneiden 
l|2001l ) for an in-depth introduction and only give a brief 
overview of the main concepts here. 

Gravitational shear 7 relates to reduced shear g as 



7 = 9- 



(10) 



where k is the projected surface mass density in units of the 
critical surface density, n = E/Ec, with 

E =^_^ an 

The latter contains a geometric factor composed of angu- 
lar diameter distances Ds,d,ds from observer to source, from 
observer to lens and from lens to source, respectively. 

The surface mass density is related to the mean tangen- 
tial component of the shear on a circle, 7t(r), by the simple 
equation 



7t(?-) = {k(< r)) -K(r) , 



(12) 



the difference between mean «: inside and on the edge of the 
circle. 

Equations [TD] and [T^] show that the observable g is in- 
variant under mass sheet transformations, k — >■ Afi:-|-(1 — A), 
a degeneracy which can be broken for instance by assuming 
a functional form of K,(r). 

In the following we use two common density profiles, 
the singular isotherm al sphere (SIS) and the NFW profile 
IjNavarro et al.lll996l ). Their densities are given by 



psis{r) = 



27rGr2 ' 

with a constant velocity dispersion a.^, and 

Po 



PNFw(r) 



(r/rs)-{l+r/r,Y 



(13) 



(14) 
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Figure 10. Tangential alignment (green circles), B mode (red 
squares) and singular isothermal sphere fit (dashed line) inside the 
interval delimited by dotted vertical lines. Note that the B mode 
is consistent with as expected for a bias-free shape catalogue. 



The NFW profile has two free parameters which are com- 
monly expressed in terms of the mass and concentration 
CAm = rAm/rs OT C/\c = T Ac / T a instead of the central den- 
sity po and the scale radius Ts used in the equation above. 

The projected surface mass density and tangential shear 
profiles of a SIS can be readily integrated as 



K.{r) = 7t(r) = 27r 



(t) 



^ Ddsl 
Ds r 



(15) 



For the NFW profile, wc refer the reader to the calculations 
in .Wright fc Braincrd. (,200a '). 



6.2 Tangential Alignment 

We first measure the tangential alignment signal around the 
brightest cluster galaxy (cf. Fig. I10|) . For this we calculate 
the mean reduced tangential shear gt in radial bins and con- 
vert it into tangential gravitational shear 7t assuming a sin- 
gular isothermal sphere (SIS) profile, for which (cf. Egns. fTOl 
andlTSll 



It = 



gt 



i + pt 



(16) 



The effective distance ratio Dds/Dg is calculated from the 
mean individual galaxy values in each bin. The SIS fit yields 
a velocity dispersion of a« = 1216 ± 76 km s~^, measuring 
in the radial range of 400 to 3000 kpc projected distance. 
As it follows from Eqn. ll3l 



r200c = 



T200r, 



V2o-„ 



(2025 ± 126) kpc and 



(17) 



v2f"u 



lOHo^n™- (1 + 2)3 



(2979 ± 184) kpc , (18) 
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which impUes a mass of M2ooc,sis = (14.0 ± 2.6) x lO" M© 
and Maoom.sis = (21.0 ± 3.8) x lO" M@. 

6.3 Significance Map 

For a first two-dimensional view of tlie lensing signal, we 
measure the surface mass density or, equivalently, tangen- 
tial gravitational shear insi de circular weig hted apertures, 
so-called aperture masses (|Schneideii 119961 ). We show the 
significance of aperture masses above zero as a function of 
position in Figure 1111 For this we use a Gaussian weight 
function 



w{\0\] 



exp(- 




my{2ai 



\e\) < 3cr„ 
otherwise 



(19) 



to calculate the significance, defined as the ratio be- 



tween aperture mass an d its 

(Bartclmann & SchncideB l200ll . 
fSchirmer et al.ll2004l ). with 



uncertainty, 
their Sec. 



5.3 and 



"■A'/a 



i^w2(|0_0^|).|^^, 



(20) 



where gi^t is the tangential reduced shear of galaxy i mea- 
sured with respect to 0. For the width of the aperture we 
use a-ui = 3arcmin. We detect the cluster as a peak in the 
aperture significance map with a significance of 5a, cen- 
tred on the BCG with a deviation of only few arcseconds. 
The significance map indicates an anisotropic distribution of 
mass around its centre, with additional peaks towards the 
north-eastern, south-western and southern direction from 
the BCG. We also check the significance of aperture B-mode 
peaks and find them to be consistent with a random field. 

6.4 Mass Density Map 

We calculate a den sity map using the method of 
iKaiser fc Squired ([1993) • For this we use all sources that sat- 
isfy the background cut of Eqn.[T] These sources have a mean 
ratio of angular diameter distances of (Dds/Ds) = 0.59. The 
K map has a pixel size of 1 arcmin and is smoothed with a 
Gaussian of cr = 1.2 arcmin width. The map is shown in 
Figure [T^] for the full field of view. 

6.5 NEW Model 

We further perfo rm a likelihood analysis (see 



ISchneider et al. 



( Navarro et al. 



2000) for a two-parametric NFW profile 
19961 ) of the halo, fitting the concentration 
parameter C200m and the mass M^m-m simultaneously. We 
use background galaxies above a projected separation of 
1 Mpc at the cluster redshift, since inside this limit some of 
the tested models predict shear above the weak regime and, 
in addition, astrophysical effects near th e core of the cluster 
disturb the expected signal most (cf. iMandelbaum et al.l 
I2OIOI ). This justifies the assumption of a model-inde pendent 
disper sion of intrinsic shapes in Eqn. 15 of Schne ider et al.l 
(|200Cl ). We use g^ ~ 0.3, the dispersion of our best-fit 
subtracted shape residuals, and determ ine confidence limits 
using the Ax^ statistics of lAvnil ([197a). 
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Figure 11. Aperture mass significance map calculated according 
to Eqns. ri9l and l20l Contour lines show the significance of aperture 
mass measured on background galaxies according to the cut of 
Eqn. [1] which is centred on the BCG at a significance of 5(T. 
The background colour shows the aperture mass significance of 
the residual shape catalogue after subtracting the best fit NFW 
model of the central halo (cf . Section I6.5| l . Towards the lower left 
is the south-eastern corner of the image. 
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Figure 12. Density map, overlayed on colour image of RXC 
J2248. 7-4431. 
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mass definition 


our NFW fit 




literature value 


source 




method 


M200m 


33.ii«:« 


28.2 iatat 3.6 ±sys 9.3 


Williamson et al. (2011") 




SZ 


Mioic 


32.2l^ej 




19 ±6 


this work. Section 14.31 and Wen et al. 


2010 


luminosity 


M2OOC 


22.815;? 




42ir 






kinematics 


M50OC 


12 7+3-'' 
^^^■'-2.6 


11 


12.25 ±0.21 

5 istat 2.6 ±sy. 0.5 


Planck Collaboration et al. (201 lb~) 
Planck Collaboration et al. f2011b') 




X-ray / XMM-Ncwton 
SZ 


M250OC 


9 0+0.8 
2-8-0.6 




5.3 ±2.6 


Comis et al. C2011) 




X-ray / Chandra 



Table 1. Best fit mass and confidence region of our NFW model (M200m) converted to other definitions of mass inside spherical 
overdensity using the best fit concentration of C200m = 2.6. All values are given in units of 10^'^Mq. 
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Figure 13. Likelihood contours of the two-parametric NFW pro- 
file of RXC J2248. 7-4431. Black, solid (blue, dotted) contours 
show the combined (projected) confidence regions (intervals) for 
M200m and C200m with no concentration prior. Red, solid (ma- 
genta, dotted) contours for the combined (project ed) confidence 
region s (intervals) using t he concentration prior of lBuUock et alj 
ll200ll ): lDuffv et al.l | |2008| ). All contours are drawn at the la and 
2(7 confidence levels. The black cross (red triangular symbol) in- 
dicate the best- fit solution without (with) concentration prior. 



1.5 
1.0' 



The la confidence limits for mass and concentration in- 
dividually are M2oom = 33.ll^;^ • IO^Mq, C2oom = 2.6 
where because of a degeneracy the very high masses only 
occur at unlikely low concentrations. Confidence contours 
are shown in Figure [13] 

Since there is prior knowledge about the concentration 
of dark matter haloes at given mass and redshift, we can 
reduce the uncertainty of our measurement by multiplying 
the likelihood with a concentration term 



p(c|M,2) = e-('°s=-'°^^(''^'^»'/(2''' 



c) 



(21) 



where we adopt th e mass-concentration relation c{M, z) of 
iDuffv et al.l (|2008[ ) and a lognormal dist ribution of concen- 
trations according to lBullock et al.l (|200ll ) with criogc = 0.18. 
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Figure 14. Projected excess mass inside a circle around the BCG 
of RXC J2248. 7-4431, shown for the best-fit NFW model (red, 
solid line) and directly measured from the observed tangential 
shear profile (black, dashed line), fixed at r200m (vertical line). 







This leads to confidence Hmits of Af200m = 30.21^;; • W^M^ 
and C2oom = 3.4lJ;g. 

For reference and comparison, we give correspond- 
ing masses in different definitions of spherical overdensity 
in Table [T] Our result is cons istent inside the error lim- 
its with the X-ray mass of the Planck Collaboration et al.l 
(l2011bl') and the SZ masses of IWilliamson et al.l l|201ll ) and 
IPlanck Collaboration et al.l(l2011bl ). The X-ray mass at high 
overdensity of I Comis et al.l l|2011l ) is higher than our best fit, 
yet relatively uncertain and likely dominated by highly con- 
centrated baryonic matter in and around the BCG, which 
is not correctly modelled by a global NFW halo. The mass 
calculated from richness and luminosity in Section 14.31 is 
significantly lower than our weak lensing result, yet likely 
contains systematic uncertainties because it was calibrated 
with different data. Spectroscopic velocity dispersion of clus- 
ter members yield s a significantly higher mass estimate 
(JGomez et al.ll2012l ). discussed in more detail in Section (T) 
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Figure 15. Left panel: Density of galaxies with \z — z^i\ < 0.06 • (1 + z^i) in our photometric redshift catalogue centred on the BCG 
position, measured as a kernel density with the minimum-variance (Epanechnikov) kernel of 1 arcmin radius. Right panel: Positions of 
red (red), moderately star-forming Sb/c-like (blue) and star-forming (cyan) galaxies around the centre of RXC J2248. 7-4431 in the same 
redshift slice. The areas delineated by black lines indicate the regions used for the SED analysis in Table[2] The two empty regions visible 
at radii of ssS arcmin from the center correspond to positions of bright stars and their associated masks. 



In addition, we plot mass enclosed in a cylinder of varying 
radius, centred on the BCG, in Figure [HI 

Deviations from the spherically symmetrical NFW pro- 
file are known to exist in several different forms, such as 
correlated secondary h aloes, asphericity of th e main halo or 
filamentary structures (|Dietrich et al.ll2012br ) and can infiu- 
ence th e accuracy of wea k lensing cluster analyses signifi- 
cantly IjGruen et al.ll2011f ). In order to investigate these ef- 
fects, we plot the aperture mass significance map on a shape 
catalogue with the signal of the best-fit model subtracted 
in Figure 1111 While the central peak is described well by 
the spherical profile, we find that there are additional peaks 
towards the north-eastern, south-western and southern di- 
rection from the BCG, with the northern one potentially 
highly elongated towards the north-north-east. This is in 
line with the density map (cf. Section I6.4|l and the shape 
of the galaxy concentration in three dimensional space (cf. 
Fig. [4]). These structures therefore most likely represent cor- 
related structures physically close to the main halo of RXC 
J2248. 7-4431. 



ON THE HYPOTHESIS OF A RECENT 
MERGER IN RXC J2248. 7-4431 



peak 


^gal 


red [%] 


mod. SF [%] 


SF [%] 


central 1' 


78 


41 ±6 


41 ±6 


18 ±4 


ENE 


67 


33 ±6 


37 ±6 


30 ±6 


SW 


50 


32 ±7 


36 ±7 


32 ±7 


O 


239 


22 ±3 


36 ±3 


42 ±3 


central 0.5' 


29 


48 ±9 


41 ±9 


11±6 


NE 


29 


51 ±9 


28 ±8 


21 ±8 



I Gomez et al.l l|2012l ) present evidence for a recent merger of 
RXC J2248. 7-4431 based on optical imaging, spectroscopy 
of 51 cluster members and analysis of the X-ray emission of 



Table 2. Fractions of galaxies by SED type for the galaxy density 
peaks identified in Section l7.1l and shown in Table |2] 



the system. We briefly recollect their findings and compare 
them with our results from deeper images with photometric 
redshift information and weak lensing. 



7.1 Galaxy number density 

The cluster member number density, selected either based 
on spectrosco pic redsh i ft or in colour-magnitude space, is 
found by Gomez et al.l (120121 ) to be multimodal (their Fig- 
ures 7-9) . We study the projected density of galaxies in the 
region of the cluster in Figure fT S] The multimodality is con- 
firmed, where the two peaks of I Gomez et al.l (|2012l ) form a 
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highly enlongated central region of approximately 2 arcmin 
diameter. An additional peak towards the east-northeastern 
(ENE) direction (no1: insid e the field of view of the images 
used by [Gomez e t al.ll2012l ) and a less dense but still visible 
peak to wards the southweste rn (SW) direction (also iden- 
tified by iGomez et alj (|2012l ) , but below their significance 
threshold). 

Our mean photometric redshifts of galaxies inside a 
cylinder of 1 arcmin radius around all three peaks are mu- 
tually consistent with each other and with the spectroscopic 
mean cluster redshift. The uncertainty of the ensemble mean 
of 5z ~ 0.005 for each peak, however, would only allow a sig- 
nificant detection of peculiar motions of several 1000 km s^^ 
along the line of sight. 

Comparing the galaxy density map with Figure [TTl we 
find a rough correspondence of the two galaxy density sub- 
peaks with the two residual peaks in the aperture mass sig- 
nificance after subtracting the best-fit model for the central 
halo. The third significant additional aperture significance 
peak towards the south does not correspond to a structure 
clearly visible in the galaxy number density. 



7.2 Galaxy SEDs 

We examine the distributions of best-fitting SED types from 
the photom etric redshift code f or the galaxy density peaks 
identified in lGomez et al.l l|2012r ) and in the previous section. 
To this end, we select gala^xies inside cylindrical volumes of 
1 arcmin (0.5 arcmin) radius and \z — Zci\ < 0.06 • (1 -I- Zd) 
around the three peaks (and the two central peaks of 
iGomez et al.ll2012l ). Results are shown in Figure [TS] and Ta- 
ble [2l discriminating red, moderately star-forming Sb/c-like 
and star-forming SEDs. 

The galaxy populations in the ENE and SW peak con- 
tain a smaller fraction of red galaxies than in the very centre 
of the cluster, yet a larger one than elsewhere at similar sep- 
aration from the BCG (in the region labelled O). 

Even if their distance from the centre in three dimen- 
sions is higher than in projection, the latter observation can- 
not be explained by the environmental influence of the core 
of RXC J2248. 7-4431 alone and indicates that ENE and SW 
are evolved neighbouring structures in the outskirts of a ma- 
jor cluster. 

For the two peaks identified bv lGomez et al.l (|2012l ). the 
galaxy populations are too small to make significant state- 
ments about differences in SED distribution (cf. Table [21 
central 0.5 arcmin and NE). 



7.3 Centroid offsets 



[Gomez et al.l ([2013) note an offset between the BCG and 
the centroid of X-ray emission when measuring the latter on 
larger radii. We note that the point of highest signiflcance 
of the lensing signal (cf. Figllip is less than 20 arcsec off the 
BCG, which is consi stent with th e random offsets expected 
from shape noise (cf. lDietrich et al. 2012a). The strong lens- 
ing model of RXC J2248.7-4431 ([Monna et al.l |2013| ) also 
shows no offset between the central galaxy and the peak of 
the projected density field. 



7.4 Discrepancy between lensing and other mass 
estimates 

Weak lensing mass estimates are not inffuenced by the as- 
trophysical state of a system, unlike estimates based on the 
dynamic state of cluster members or the intra-cluster gas. 
Therefore a discrepancy between weak lensing and other 
mass estimates could suggest the infl uence of dynamic as - 
trophysical processes such as mergers. I Gomez et al.l (|2012[ ) 
hypothesise that a lensing mass signiflcantly lower than the 
dynamical and X-ray masses found for RXC J2248. 7-4431 
could help confirm merger ac tivity. 

The mass calculated by [Gomez et al.l (|2012[ ) based on 
the velocity dispersion of cluster members is slightly in- 
consistent with our weak lensing mass on the high side 
(M2oVc = 42li^TO"Af0 versus M^tc = 22.8ll^lO"M0), 
indicating that opposite bulk motion of different galaxy pop- 
ulations along the line of sight could be present. Whether 
there is indeed substruct ure in the velo city fields of the clus- 
ter members, for which IComez et al.l (2012) find marginal 
evidence will be investigated in detail with the forthcom- 
ing VLT-CLASH large spectroscopic programme (P. Rosati, 
private communication) . 

We see no significant discrepancy, however, of weak 
lensing with X-ray or SZ mass estimates of the cluster. 
I Gomez et al.l (|2012l ) find an inconsistency between masses 
inside 110 kpc apertures from X-ray modelling in hydrostatic 
equilibrium and a strong lensing model. The latter, how- 
ever, is based on a candidate multiply imaged system iden- 
tified in their relatively shallow photometry which deeper 
observations have r evealed to be not from the same source 
(JMonna et al.ll2013r ). At this point, there is therefore also no 
compelling evidence of a discrepancy between lensing and 
X-ray mass estimates in the core of RXC J2248. 7-4431. 



7.5 Influence of neighbouring structures on 
lensing mass 

We attempt to quantify the influence of the structures in 
the environment of RXC J2248. 7-4431 on its weak lensing 
mass estimate by fitting a model with multiple haloes to the 
observed shear signal. To this end, we place two halo es with 
fixed concentration-mass relation ([Duffv et al.ll2008i ) at the 
position of the BCG and the centre of the east-northeastern 
peak (cf. Section [7T). The confidence region for the mass of 
the secondary peak is in this case M§)^^ = 4.2^ jg • 10^*Mq . 
Subtracting the best- fitting model of the secondary peak and 
fitting the central halo of RXC J2248. 7-4431 with two free 
parameters yields M2oom = 24.0lg'g ■ 10"'^'' Mq and C2oom = 
2.8^12. The corresponding mass for comparison with other 
estimates is Msooc = Q-Gtl.l ' IO^Mq. 

We note that this lower mass of the central peak is still 
in agreement with the SZ and X-ray masses listed in Ta- 
ble [l] Is is unclear whether or not the mass of the secondary 
peak should in fact be included in a M2oom estimate of RXC 
J2248. 7-4431 since the projected separation of 3.5 arcmin is 
well below r2oom of the central system. While this exem- 
plifies the sensitivity of clust er weak lensing to correlated 
structures (|Gruen et al.ll201H ). it does not change the con- 
clusions of Section 17.41 
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8 SECONDARY CLUSTER AT Z ^ 0.6 

In redshift space density maps we detect a secondary peak 
around {a,S) = (22''49'"37.1% -44°43'04"), at the position 
of a bright early-type galaxy with photometric redshift of 
z = 0.66 and corresponding absolute rest-frame magnitude 
of Afij^voga ~ —25.1 which appears to be at the centre 
of another cluster of galaxies. At a separation of 14.5 ar- 
cmin from the centre of RXC J2248. 7-4431, this configura- 
tion of lenses is significantly more se parated than the c ases 
of clustcr-cluster-lensing discussed in lZitrin et al.l ([2013). In 
the following, we present a weak lensing and a simple strong 
lensing analysis of the newly discovered cluster. 

The mean redshift of the 12 closest, visually colour- 
selected member galaxies is z = 0.58, which is the redshift 
we use for the following weak lensing analysis. After sub- 
tracting the best-fit NFW model of the central peak of RXC 
J2248. 7-4431 (cf. Section[S3I, we model this structure as an 
NFW halo positioned on the BCG using the concentration 
prior of Ean. l21l and a radial range of 0.5 Mpc ^ r ^ 3 Mpc. 
The confidence contours are shown in Figure [TB] As a re- 
sult of the low background number density and location of 
the cluster near the edge of the image, the mass is only 
weakly constrained and not detected at the 2cr level. We 
find the confidence region to be M2oom = 4.0^26 ' IO^^Mq 
at a concentration of C200m = 4.31^4. Owing to the lack of 
completeness at the required depth, a richness and luminos- 
ity derived mass estimation (cf. Section r4.3|) is not feasible. 
Figure \T7\ shows a colour image of the central part of the 
system. 

We identify four blue sources in a cross-like orienta- 
tion around the BCG with separations of ~ 8 arcsec and 
~ 17 arcsec along the small and large axis, respectively, as a 
candidate of multiple imaging. Assuming a cluster redshift 
of zi = 0.58, the redshift of the blue source at its photo- 
metric value of Zs = 1.58, and a halo with projected ellip- 
ticity and orientation following the light distribution of the 
BCG, we can build a strong lensing model. For a nonsingu- 
lar halo with 4 arcsec (27 kpc) core radius, we get a mass of 
~ 1.4 • IO^^Mq within a radius of 6.8 arcsec, the Einstein 
radius at this redshift. 



9 CONCLUSIONS 

We present the first weak lensing analysis of the cluster of 
galaxies RXC J2248. 7-4431. We calculate photometric red- 
shifts of background and cluster galaxies and determine the 
richness and total luminosity of the system. We find the mass 
and concentration of RXC J2248. 7-4431 in a NFW hkelihood 
analysis to be equal to M2oom ~ 33.1^.^'^ ■ 10^*Mq, in good 
agreement with previous X-ray and SZ mass estimates, and 
C200m = 2.6l]^Q, which is in the lower range of expected con- 
centration for a system of the given mass and redshift. The 
subtraction of the best-fitting model leaves three marginally 
significant aperture mass peaks in the vicinity of the main 
halo, which likely correspond to substructure and surround- 
ing structures of RXC J2248. 7-4431. We co nfirm and add to 
some of the evidence given bv lComez et al.: ( 2012) for an on- 
going major merger of the system in terms of the multimodal 
galaxy density field and a discrepancy between (low) weak 
lensing and (high) dynamical mass estimate. We do not, 
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Figure 16. Likelihood contours of the two-parametric NFW pro- 
file of secondary cluster with z = 0.58. Red, solid (magenta, 
dotted) contours for the combined (projected) confidence regions 
(i ntervals) for M2oom a nd C20Qm using the concentration prior 
of iBulIock et alJ l|200lh : IPuffv et all l|2008h . Both contours are 
drawn at the Ic and 2cr confidence levels. The triangular symbol 
indicates the best-fit solution. 




Figure 17. Colour image of second cluster at {a, 5) = 
(22'^49™37.1=, -44°43'04"), z fs 0.6, with a size of 2x2 arc- 
min. Note the redder colour of BCG and cluster member galax- 
ies compared to RXC J2248.7-4431 (cf. Figure [U and the four 
symmetric blue sources along the major and minor axes of the 
BCG (the former ones highly blended with the BCG itself) at 
%hot ~ 1-5 • • ■ 1-8- 



© 2013 RAS, MNRAS OOO.ITHTsl 



14 D. Gruen et al. 



however, find a tension between X-ray hydrostatic masses 
and the weak lensing mass of the cluster. This remains to 
be true even when a secondary peak of galaxy density at 
3.5 arcmin separation of the core of RXC J2248. 7-4431 is 
independently included in the weak lensing model. 

We detect a second cluster at z ~ 0.6 inside our field 
of view, whose weak lensing mass is weakly constrained 
at M2oom = 4.0I2 6 ' IO^^Mq with a concentration of 



C200r7 



4.3 



-1-2.2 



In this system a strong lensing analysis 



of a candidate multiply imaged source is possible. 
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